Modulations to molecular high order harmonic generation by electron de Broglie wave 
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We present a new theory that the molecular high order harmonic generation in an intense laser 
field is determined by molecular internal symmetry and momentum distribution of the tunneling- 
ionized electron. The molecular internal symmetry determines the quantum interference form of 
the returning electron inside the molecule. The electron momentum distribution determines the 
relative interference strength of each individual electron de Broglie wave. All individual electron 
de Broglie wave interferences add together to collectively modulate the molecular high harmonic 
generation. We specifically discuss the suppression of the generation on adjacent harmonic orders 
and the dependence of molecular high harmonic generation on laser intensities and molecular axis 
Q_< alignment. Our theoretical results are in good consistency with the experimental observations. 

PACS numbers: 42.65.Ky, 32.80.Rm 

Laser-molecule interaction has become an active research area in high field physics. The molecular multi-center 
structure strongly affects its interaction with the external laser field, resulting in abundant new physics [11, H, H, 0, 
JlG, EH, EH EH ■ Recent studies show that the molecular orbital symmetry strongly affects the single 
ionization [lj, |2|, |3|; the nuclei motion in molecular dissociation and Coulomb explosion can be monitored by studying 
the electron rescattering process 0, the molecular orbital can be imaged by studying the angular distribution of 
coincident ion pairs in the molecular double ionization [6(; by studying the spectral profile of molecular high order 
harmonic emission, one is able to reconstruct the tomography of the molecular structure 0, EH EH EH • Of particular 
interest is the molecular high order harmonic generation (HHG). The molecular multi-center structure provides the 
possibility to actively modulate the molecular HHG to generate a single attosecond pulse. 

The ionization process of a molecule in an intense laser field is described using the well-known simple-man's 
model[l4j]. First one molecular electron tunnels through the potential barrier assisted by the laser field, then the 
, electron in the continuum state may return to the molecular ion as the laser field reverses its phase. Different from 
■^j- ■ the rescattering process in the laser-atom interaction, the multiple atomic sites inside the molecular ion provide mul- 
tiple options for the electron recombination. The path-differences of the returning electron to the multiple atomic 
sites result in the interference of electron de Broglie wave (e lk r ) which significantly affects the molecular ionization 
and HHG. This was first noticed by Lein and his colleagues in studying the interaction of laser with molecules . 
Lein et al. showed that when recaptured into the bonding orbit of a molecule, the electron waves from the two atomic 
site emitters (with positions ri and ra) interfere as, 

e lkri +e lk - r2 , (1) 
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which results in constructive or destructive interference when k • R (R = ri — ra) is a multiple or a multiple and a half 
of the laser wavelength. As a result, the n th -ordei molecular harmonic emission (nui = Ek, Ek — k 2 /2 ) is enhanced 
(or suppressed). 

Recently two experiments were carried out independently to study the HHG of CO2 in an intense laser field with 
similar experimental configurations only differing in laser intensities. The suppression of HHG was observed for a 
range of adjacent harmonic orders in both experiments [Hil Il7j|: while the suppressed harmonic spectral ranges were 
displaced by about 8 harmonic orders with respect to each other. This cannot be explained by Eq. |T]). 

In this letter, we present a new theory based on the modified molecular Lewenstein model to describe the laser- 
molecular interaction. We show that similar to the Young's two-slit interference, the returning electron (tunncling- 
ionized in the external laser field) wave interferes on the multiple atomic sites inside the molecule ion. The returning 
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electron generally has a momentum (k) distribution. Each single electron momentum state (the electron de Broglie 
wave) self-interferes on the multiple atomic sites. All interferences of different electron momentum states (k) add to- 
gether to affect the molecular ionization and HHG. The interference structure is uniquely determined by the molecular 
structure. The electron momentum distribution is a function of laser intensity and molecular axis alignment angle in 
the external laser field. Thus the molecular HHG changes as laser intensity or molecular axis alignment varies. We also 
show that molecular HHG at adjacent harmonic orders has similar dependence on electron momentum distribution, 
thus the suppression or enhancement on molecular HHG occurs over a broad range of adjacent harmonic orders. Our 
theoretical results agree well with experimental observations. 

To directly compare with experimental observations, we study the HHG of CO2 which has anti-bonding orbital as 
the highest-occupied state and can be approximately considered as a two-center system when it interacts with the 
external laser field [l6l. [l7|. Using the molecular Lewenstein model, the dipole acceleration of CO2 to yield the HHG 
radiation in an external laser field (polarized along x) is written as[l8l [l9l I20I l2l| . 

a x (t) = i J d 3 k|$,(k)| 2 ^ dt'A(t') 



x E(t)(2exp[-iS (k,t,t')] - exp[-iSi(k, t, i')] - exp[-iS , 2 (k, t, t')]) + c.c, (2) 

where A(t) and E(t) are vector potential and electric field of the laser pulse. Beside the action So(\t,t,t') — 
/ t *,(ft"(ik 2 -k- A(t")+I p ), two additional action terms Si(k, t, t') = 5 (k, t, t')-k-R and S 2 (k, t, t') = S (k, t, t')+k-R 
are introduced accounting for the effects induced by the two O-atomic sites. The three phase terms carrying these 
actions interfere with each other and modulate the molecular HHG(l9l. l20l . l2ll [22J. ^(k) is the Fourier transfor- 
mation of the atomic wavefunction tfi under linear combination of atomic orbitals-molecular orbitals (LCAO-MO) 
approximation in which the molecular wavefunction is expressed as tpi — tpi(r, — R/2) — tfi{v, R/2). 

After adopting pole approximation [23] and performing time- integral in Eq. ([2]), the Fourier transformation of Eq. 
© gives the HHG spectral amplitude of (702 0,121 

S(n) oc J2 K^''* " e-^-%)MV)\ 2 Ji(^^)J m (^), (3) 
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where J;() is the I th order Bessel function, m and I are numbers of photons that the electron absorbs or emits at 
ionization or recombination, with n = m±l + l. The electron's momentum is now given as fc 2 /2 = (m±l)o;— I p , where 
I p is the ionization potential of the molecule. Considering the acceleration effect by bound potential when the electron 
is in the vicinity of the core, the modified electron momentum k' (parallel to fc) is given as fc' 2 /2 = fc 2 /2 + Z^j^. 

Eq. §5§ shows that the molecular HHG is completely determined by the structural factor( e lk ' "a — e~' k ' which 
corresponds to the molecular internal symmetry (here it describes the anti-bonding orbital of CO2) and the momentum 
distribution [the rest part of Eq. [[5])] of the rescattering electron which is a function of the laser intensity and molecular 
axis alignment angle (with respect to the direction of the external laser field). Eq. ([3]) also shows that the electron 
momentum distribution and consequently the molecular HHG spectral profile vary for different laser intensity and 
molecular alignment angle. In the following, we numerically study the HHG of CO2 under various experimental 
conditions and compare the results with available experimental data. 

We first study the contribution to molecular HHG from electron momentum states with momentum k parallel to 
E, because they were considered to contribute the most to the HHG. As shown in Fig. 1, for both harmonic orders 
S(n — 13) [Fig. 1(a)] and S(n — 23) [Fig. 1(b)], contributions from individual electron momentum states (with 
k parallel to E) are significantly suppressed for smaller k but remain unchanged for big k, when considering the 

interference effect that is due to the anti-bonding symmetry of CC>2- An integral of J G?f2 J Q k 2 dk was then conducted 
to consider contributions from all electron momentum states. Both S(n = 13) [Fig. 1(c)] and S(n — 23) [Fig. 1(d)] 
are constant for big k, indicating that the contribution from big-fc electron momentum states is negligible. When 
considering the interference effect, S(n = 13) is relatively suppressed while S(n = 23) is relatively enhanced. 

Next we numerically study the variation of HHG spectral profile of CO2 as a function of laser intensity and molecular 
alignment angle 9. For fixed laser intensity, the HHG spectral profiles vary for different values of 9 as shown in Figs. 
2(a) and (b) in which the plots are generated with / = 1.5 x 10 14 W/cm 2 and I = 2 x 10 14 W/cm 2 , respectively. Eq. 
([3]) shows that harmonic emissions at adjacent orders have similar dependence on electron momentum contribution. 
Consequently similar modulations occur on adjacent HHG orders. As shown in Fig. 2(c), with I = 2 x 10 14 W/cm 2 
and 9 — 30°, HHG is suppressed at a range of adjacent harmonic orders centering at n = 25. The center shifts to 
n = 35 for I = 3 x 10 14 W/cm 2 [Fig. 2(d)], qualitatively agreeing with recent experimental observations (l6l [l7l ]. 
(Intensity stabilization and calibration of intense ultrashort laser pulse remain to be a challenge in the experiment.) 
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To clearly exhibit the dependence of molecular HHG on the laser intensity, we numerically study the HHG at 
two harmonic orders. For n — 13, our numerical results show that the harmonic generation is suppressed at I = 
1.5 x 10 14 W/cm 2 [Fig. 3(a)] and enhanced at I = 2 x 10 14 W/cm 2 [Fig. 3(b)]. For n = 23, we show the suppression 
of harmonic generation at / = 2 x 10 14 W/cm 2 [Fig. 3(c)] and enhancement at / = 3 x 10 14 W/cm 2 [Fig. 3(d)]. 

Lastly, it is interesting to find out under what kind of condition the single electron momentum picture given by 
Lein et al. may be used in predicting the suppression position in the HHG spectrum. Examining Eq. (|2|) shows 
that if k • R (describing the molecular structure) is much smaller than the action Sq (corresponding to the electron 
trajectory), then Sq can be extracted out as a common factor. 

The Fourier-transformation of Eq. ([2]) can now be simplified into Eq. {T]). Assuming A(i) = ^-cos^t), we have 

- E 2 

Sq ~ Ek.T, where T is the period of the laser field and E^is the average kinetic energy of the electron with Ek ~ 2 4 ^ 2 
and k ~ Eq/co. Then we reach such a relation 

R < a = E q /lo 2 , (4) 

where a is the electron's quiver amplitude in the laser field. Eq. ((4| shows that if the electron's quiver amplitude is 
much larger than the molecular internuclear distance, or in other words, the time difference for the rescattering electron 
to recombine with the different atomic sites are negligible comparing to the time for electron being in the continuum 
state, then Eq. ([I]) can be used to approximately describe the molecular HHG process. Comparing the internuclear 
distances R = 2.0 a.u. for and R = 4.4 a.u. for CO2 with the electron's quiver amplitude of a = 21 a.u. in a laser 
field with intensity of / = 1.5 x 10 14 W/cm 2 , it is easy to understand that Eq. (JTJ) gives a better estimation to the 
suppression of HHG of -ff 2 + [Ii than to the suppression of HHG of COi. 

In conclusion, we develop a new theory about the molecular HHG in an intense laser field. We illustrate that the 
molecular HHG at each order is contributed by a group of interference effects on the Young's slit (multiple-atomic 
sites) with each produced by an independent electron de Broglie wave (eigen electron momentum state) and weighed 
by the electron momentum state amplitude. We show that the single electron momentum picture illustrated by Lein 
et al. is only an approximation of our theory at high laser intensity where the electron's quiver motion amplitude 
is much larger than the molecular internuclear distance. Our theoretical analysis shows that the difference in the 
observed suppressed HHG spectral ranges in recent experiments is attributed to the different laser intensities applied 
in the measurement. We further suggest that the molecular HHG in an intense laser field may be generally described 
in the form of Eq. ([3]) by simply replacing the interference factor with the appropriate structure factor of the molecule 
in use. 
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Figure 1: Calculated distributions of contributions to the HHG of CO2 at n = 13 (a) and n = 23 (b) from individual electron 
momentum state with k parallel to E, considering (black open dots) and without considering (red filled dots) interference 
effect. Integrated contributions from various electron momentum states (J dfl J k 2 dk) to the HHG of CO2 at n = 13 (c) and 
n — 23 (d), considering (black open dots) and without considering (red filled dots) interference effect. I — 1.5 x 10 14 W/cm 2 
and A = 800nm. (The double-line-structure is due to the double values of k with k 2 /2 = (m ± 1)lu — I p .) 
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Figure 2: Calculated HHG power spectra of CO2 at different molecular alignment angles with (a) / = 1.5 x 10 14 W/cm 2 
and (b) I = 2 x 10 14 W/cm 2 . Calculated HHG suppression (spectra normalized to that averaging over all angles) with (c) 
/ = 2 x 10 14 W/cm 2 and (d) I = 3 X 10 14 W/cm 2 . A = 800nm. 
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Figure 3: Integrated contributions from various electron momentum states (f dQ. k 2 dk) to the HHG of CO2, considering 
(black open dots) and without considering (red filled dots) interference effect, (a) n — 13 with I = 1.5 x 10 14 W/cm 2 and (b) 
n = 13 and I = 2 x 10 14 W/cm 2 , (c) n = 23 with 7 = 2 X 10 14 W/cm 2 and (d) n = 23 and I = 3 x 10 14 W/cm 2 . A = 800nm. 



